130 GeV dark matter and the Fermi gamma-ray line 
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Based on tentative evidence for a peak in the Fermi gamma-ray spectrum originating from near 
the center of the galaxy, it has been suggested that dark matter of mass ~ 130 GeV is annihilating 
directly into photons with a cross section ~ 24 times smaller than that needed for the thermal 
relic density. We propose a simple particle physics model in which the DM is a scalar X, with a 
coupling Axels'! 2 to a scalar multiplet S carrying electric charge, which allows for XX — > 77 at 
one loop due to the virtual S. We predict a second monochromatic photon peak at 114 GeV due 
to XX — ► yZ. The S is colored under a hidden sector SU(N) or QCD to help boost the XX —¥ 77 
cross section. The analogous coupling Xhh 2 \S\ 2 to the Higgs boson can naturally increase the partial 
width for h — > 77 by an amount comparable to its standard model value, as suggested by recent 
measurements from CMS. Due to the hidden sector SU(N) (or QCD), 5* binds to its antiparticle 
to form S-mesons, which will be pair-produced in colliders and then decay predominantly to XX, 
hh, or to glueballs of the SU(N) which subsequently decay to photons. The cross section for X on 
nucleons is close to the XenonlOO upper limit, suggesting that it should be discovered soon by direct 
detection. 
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Refs. [1, 2] have recently found tentative evidence for 
a narrow spectral feature at E 1 = 130 GeV in the Fermi- 
LAT [3] data (a 4.6cr excess, or 3.3cr taking into account 
the look-elsewhere effect), and have interpreted it as pho- 
tons from the annihilation of dark matter (DM) of the 
same mass. The Fermi collaboration does not yet re- 
port such a signal, but their most recent upper limit of 
(av) ~ 10 _27 cm 3 s _1 (assuming an Einasto profile) for 
130 GeV DM to annihilate into two photons [4] is con- 
sistent with the required cross section found in [2]. The 
DM interpretation was bolstered in rcf. [5] , which showed 
that the two-photon annihilation channel gives a better 
fit to the feature than do other final states leading to 
photons, the others tending to give a broader peak than 
is observed. Ref. [6] has suggested that the excess has 
an astrophysical origin associated with the Fermi bub- 
ble regions, but rcf. [5] claims to locate the spatial re- 
gions in which the signal is maximized, indicating that 
the strongest emission is coming from close to the galac- 
tic center and not the Fermi bubble regions. In this note 
we adopt the annihilating DM hypothesis and propose a 
model which can account for the monochromatic photon 
line. 1 





FIG. 1: Feynman diagrams for the annihilation XX 
mediated by virtual S. 



TO 



'Electronic address: jclme@physics.mcgill.ca 

1 For an alternative model involving an extra U(l) gauge boson 
see [7]. See also [8] for an earlier model that can provide gamma 
ray lines from DM annihilation. 



Dark matter (here denoted by X) should couple only 
weakly to photons, if at all, at tree-level [9, 10]. One way 
to insure the "darkness" of the DM is for it to couple 
to photons only via loops. At one loop, the DM should 
couple directly to charged particles S. To make a renor- 
malizablc coupling of this type, both X and S must be 
bosons, since the stability of X and the conservation of 
charge require X 2 and IS"! 2 . This leads us to consider the 
interactions 



X 2 \S\ 2 + \ h \H\ 2 \S\* 



m 2 x 2 



(i) 



between X, the Higgs doublet H, and S. The second cou- 
pling is not necessary, but neither is there is any reason to 
forbid it, and in fact we will show that it can naturally 
give rise to an interesting enhancement in the h — > 77 
branching ratio, for the same values of the S mass and 
charge as needed to explain the Fermi line. The third 
coupling is useful for achieving the correct relic density 
of X [11], as we will discuss. The stability of X is insured 
by the Z2 symmetry X — > — X. 

Decays of S. It is necessary to make S unstable 
in order to avoid charged relics, on whose abundance 
there are very stringent bounds from terrestrial searches 
for anomalous heavy isotopes [12, 13] and from their ef- 
fects on big bang nucleosynthesis [14, 15]. We will also 
find it useful to let S transform under QCD or a hid- 
den SU(N) gauge symmetry, in order to boost the cross 
section for XX 77. Suppose S is in the fundamen- 
tal representation of SU(N) for definitencss. If SU(N) 
is QCD and S has charge 4/3, it can decay into right- 
handed up-type quarks through the renormalizable op- 
erator e a ^S a u R ,^u R ^.. If the SU(N) is exotic, then S 
could decay into a lighter, neutral fundamental repre- 
sentation field T and two charged right-handed fermions 
through a dimension 5 operator. For example, if S has 
charge q$ — 2, the decay into T + e + + e + occurs via the 



operator 



T* <?<V C p 

M R R 



(2) 



which could arise from a renormalizable theory by inte- 
grating out a heavy colored fermion N a carrying charge 
1. For q s = 1, there is an analogous operator to (2) in- 
volving left-handed lepton doublets, T*S a L c e a2L ei that 
mediate the decay S — > Te + D. We will focus on the 
example (2) because the larger charge qs = 2 helps to 
increase the XX —> 77 cross section, and the decays into 
electrons can lead to interesting collider signatures. Of 
course, higher-generation leptons can also appear in ad- 
dition to electrons through analogous operators, as well 
as lepton flavor-violating versions, whose contribution to 
the rare process fi — > e7 is suppressed by two loops and 
two powers of M. 

Obviously the lifetime of S depends upon the scale M 
of the heavy N a particle. Estimating the decay rate for 
S —> Tee as m|/(167rM 2 ) and demanding the lifetime to 
be less than 10 3 s [14], we find the limit M < 10 15 GcV. 
We will show that the relic neutral particle T a binds into 
stable "baryons" that make a small contribution to the 
total dark matter population. 

For simplicity we have assumed that additional cou- 
plings of T such as X XT X 2 \T\ 2 and A sr |S| 2 |T| 2 arc small. 
The former could provide a significant annihilation chan- 
nel XX -> TT* if m T < m x /2 (the factor of 1 /2 coming 
from the fact that each T must hadronize into TT* bound 
states) and if A XT is sufficiently large, while the latter has 
no particular impact on the points that follow. 

Annihilation to two photons. The model parame- 
ters relevant for the Fermi line are Xx , the mass mx , the 
charge qs (in units of e), the mass mg, and the number of 
colors N c of QCD or the hidden SU(N) gauge group. The 
annihilation cross section corresponding to the diagrams 
of fig. 1 is given by 



E\M\ 2 
647rm 2 Y 



(3) 



where the squared matrix element, summed over photon 
polarizations, is 



E\ M \ 2 = ^^xN 2 c r'Al(r) (4) 

with t = m\ /m s and 

tA (t) = 1 - T^arcsh^V?) 2 (5) 

for t < 1, which we presume to be the case. Eqs. (4,5) 
can be deduced by comparing to the well-known result 
for h — > 77 from fig. 2, in which r — > m|/4m|; see for 
example ref. [16]. 

Ref. [2] determined that (av) should be approximately 
0.042 in units of the thermal relic density value (av)o = 
1 pb-c, in order to explain the Fermi gamma-ray line. 
(Version 1 of ref. [5] found a larger value, comparable to 
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FIG. 2: Feynman diagrams for the decay h — >• 77 mediated 
by virtual S. v = 246 GeV is the Higgs VEV. 




m s (GeV) 



FIG. 3: Value of S charge qs (in units of e) as a function of 
ms, needed to obtain the X —¥ 77 cross section of 0.042 times 
the standard relic density value (crv)o, assuming mx = 130 
GeV. 



(<tv)q, but this was due to an error that has now been 
corrected.) Taking this value for (av) and m x = 130 
GeV, we can find the relation between qsV^xN c and 
ms, shown in fig. 3. 

From fig. 3 we see that even if the coupling is rather 
large, ~ 3, and N c = 3, the S charge is typically 
greater than 1 (in units of e), and only reaches 1 for ms 
close to m x . 2 At the other extreme of qs — 6, the corre- 
sponding interaction strength q 2 s a = 0.26 would still be 
under perturbative control, indicating that models with 
ms up to at least 400 GeV are viable. On the other hand, 
the rate goes like qgX x , so with qs = 6 one could alter- 
natively lower A^y from 3 to 0.08 by taking ms = m X - 
Thus it is not necessary to invoke a very large value of 
Xx- 

Electroweak Precision Constraints. S must in- 
herit its electric charge from weak hypcrchargc, and 
therefore it couples also to the Z boson with strength 
qsetanOw , where 9w is the Weinberg angle. Such a 
particle, if neutral under SU(2)^, is unconstrained by 
precision electroweak constraints since its contribution 
to the p or T parameters vanishes identically [17], and it 



2 This limiting case might be of interest because it allows for 
the possibility of embedding S into an SU(2)^ doublet with the 
standard hypcrchargc assignment for extra Higgs doublets (see 
arxiv version 1 of this paper for more details). However this is 
not necessary, and we will focus on the case where S is a singlet 
of SU(2)^ and carries only weak hypcrcharge. 
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contributes only to the Y parameter [18] which is weakly 
constrained. 

Annihilation to Z bosons and photons. As noted 
above, S necessarily couples to the Z boson as well as 
to photons. Comparing the seagull vertices for |<S'| 2 77 
and |S'| 2 7Z, we can deduce that the cross section for 
XX — > Zj is related to that for XX — > 77 by the factor 



f v \ xx ^ z = 2tan 2 g w (l- 
{av)xx^ 11 y 



Am 2 x 



1/2 



= 0.56 (6) 



taking into account the reduced phase space for identical 
particles in the case of XX — > 77. Since the former 
process produces only one photon, its intensity will be 
0.28 times that of the 2-photon line. The energy of this 
single photon is given by 



E~ ( = mx — 



Am 



= 114 GeV 



x 



(7) 



We therefore predict that the spectral feature will re- 
solve into two peaks separated by 16 GeV in energy. 
The relative strength of the peaks could be modified by 
giving a different SU(2)l assignment to S so that both 
components of the doublet become electrically charged 
(see footnote 2). The current Fermi/LAT limit on 
(&v) x x ->■-{ z is 2 x 10~ cm 3 s~l [3] is compatible with 
our value. 

Relic density. The cross section for XX —> 77 
is well below that which is needed to obtain the right 
relic density, but this can still be achieved using the 
XX — > hh, WW, ZZ channels mediated by the interac- 
tions jAhx^X 2 , ^XhxvhX 2 , Xvh 3 (from the standard 
model Higgs potential) and the gauge couplings of the 
Higgs, with intermediate Higgs h in the s channel. [11]. 
The cross section for XX — > hh is 



(<rv)hh 



A hX 

64-7T m 2 x 



Th ^hx . 1 + r h /2 



2 A l-r h /4 



where 77, = rri^/m 2 ^ = 0.94 and A = 0.13 assuming 
126 GeV. For the WW final state, we find 



'V, 
(8) 




(9) 

where r w = m^/m^. The cross section (<rv) zz is 
the same as {<jv)\ VW with the replacement v w — y t z = 
m %/ m x anc * a factor of 1/2 for identical particles in the 
final state. Demanding that (eru)tot = {wo) a implies that 
|A/ix| = 0.051. If there are other significant annihila- 
tion channels, this contribution must be correspondingly 
reduced so that in general \\hx\ < 0.05. For example 
if S carries QCD color, then XX — > gluons resulting in 
hadrons can be significant. There is also the possibility 
of XX -> TT* if the coupling A XT V 2 |T| 2 is sufficiently 
large and if mj- < mj/2. 

Direct detection. The ^Xhxh 2 X 2 vertex gives rise 
to the trilinear interaction XuxvhX 2 from electroweak 
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FIG. 4: Ratio of new and SM contributions to the h — > 77 
decay width as a function of ms assuming the relation in fig. 
3, m h = 126 GeV, m x = 130 GeV, \ x = 1, N c = 3, and 
A/i/Ax = 0.5, 0.75, 1 as indicated. 



symmetry breaking. (Note that the neutral component 
of the Higgs doublet is H° = -^{h + v)). The Higgs can 
therefore mediate scattering of X on nucleons V. The 
cross section for XN — > XN elastic scattering is [19] 



P Kx m N 
^ m h m X 



(10) 



where fmjsi/v is the Higgs- nucleon coupling, with / = 
E q=u ,d,sfT q + (2/9)/? G = 0.319 [20]. Using the con- 
straint X^x ^ 0.05 from the preceding relic density de- 
termination, we can evaluate (10) to find a < 1.3 X 10~ 45 
cm 2 . This is an order of magnitude lower than the 2011 
XenonlOO 90% c.l. limit of 1.2 x 10" 44 cm 2 at m x = 130 
GeV [21], but only 3 times lower than the new limit of 
3x 10~ 45 cm 2 that was recently announced [22]. Thus the 
model could be confirmed or ruled out in the near future 
by anticipated improvements [23] in the direct detection 
limit. 

Implications for Higgs decays. Because of the close 
similarity between the diagrams of figs. 1 and 2, there is 
a simple relation between {av) and the extra contribu- 
tion to h — > 77, which is especially interesting in light of 
the recent observation by CMS of an upward fluctuation 
in that branching ratio, relative to the standard model 
expectation [24], assuming of course that the indications 
of discovery of the Higgs boson with nih — 126 GeV are 
borne out [25-27]. Specifically, the squared matrix ele- 
ment for h — > 77 is related to (4) by replacing Xx — > X^v 
(where v is the Higgs VEV) and t —¥ f = m\jAm 2 s . The 
contribution of the charged scalar interferes construc- 
tively with that of the SM if A^ > [16]. To give an idea 
its relative size, we can express the extra contribution to 
the partial width of h — > 77 (here ignoring interference 
effects) in terms of its ratio to the SM contribution (see 
eq. (5), 



SM 



9^A 2 V 2 



2y/2v 2 



A 2 SM G F ; 



f 2 A 2 {f) (11) 
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where the SM amplitude is given by Asm — ™6.52 for 
to/j = 126 GeV. For given values of X^/Xx, N c and mg 
(and assuming that nix = 130 GeV) the ratio (11) is 
fixed by the Fermi-LAT cross section, and is 0(1) for a 
wide range of ms if Xh ~ 0.5Xx- We plot it as a function 
of ms in fig. 4 for several values of Xh/Xx- 

Confinement of S particles. Considering the case 
where q s = 2 and S decays arc mediated by a dimension 5 
operator (2), generically one would expect the S lifetime 
to be much larger than the hadronization time, 1/A for 
a gauge theory with confinement mass scale A. Thus 
any S particles produced in a collider will have time to 
form "mesonic" or "baryonic" bound states with other S 
or T particles. We first discuss the SS*, ST* and TT* 
mesons, which we denote by <ps, tts and 4>t respectively. 

Through the interaction (2), tts can decay into two 
electrons with a relatively long lifetime. But <f>s decays 
much faster by its constituents annihilating into XX, hh, 
77, ZZ and gg final states, where g is the gluon of the 
SU(N) theory. Since there are no lighter colored parti- 
cles in the theory, these gluons hadronizc into glucballs, 
which we denote by CI. (CI is unstable to decay into pho- 
tons, as we discuss below.) The partial decay widths can 
be estimated using T ~ n(av) where the density n is 
the square of the 5- meson wave function at the origin, 
n ~ (a'ms/2) 3 , a' is the strength of the SU(N) gauge 
interaction at the scale ms, and a is the cross section for 
SS* scattering into the desired final state. In this way 
we find the partial widths 



0S 
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64-7T 
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4>s -> 77: zz 



CI CI 



(12) 

Because of the large coupling Ax, the invisible width for 
decays into XX is typically the most important; even if 
the mass splitting mj — mx is small, for example 3 GeV, 
Txx = 170 T 77 when Aa" = 3. However, the branching 
ratio into glueballs Cl can also be large and even domi- 
nant, depending upon the unknown value of a' . The glue- 
balls can decay into photons via the Eulcr-Hcisenberg in- 
teraction FFGG/mg induced by an S loop. If the SU(N) 
confinement scale A is below the weak scale, then pho- 
tons will be the only kinematically available final states 
for the glueballs to decay into. The partial width for 
4>s — > hh can be signficant if Xh ~ Xx ■ 

The 4>t meson can decay into all the same final states 
as the (f>s by first going into two gluons of the SU(N) the- 
ory, which turn into two other particles by going through 
a loop of S. These amplitudes thus occur at two loops. 
They can be dominated by a one-loop contribution due 
to the possible interaction A|5| 2 |T| 2 . Thus we expect the 
4>t partial widths to be proportional to those of <f>g, with 



an extra loop and coupling constant suppression. 

The S and T particles can also bind together into 
baryons of the SU(N). For example if V = 3, we have 
the SSS, SST, STT and TTT baryons, which can decay 
to states with fewer S's via S —> T + 2e until reach- 
ing the stable, neutral TTT state. The former would be 
stable charged relics in the absence of this decay mech- 
anism. The latter is a dark matter candidate, but its 
annihilation cross section is too large for it to provide 
a significant contribution to the total DM density. The 
TTT annihilation cross section is of order a 12 /(imr) 2 , 
which for a' = 0.1 (similar to the strength of QCD at 
these energy scales) is some 2 orders of magnitude larger 
than needed to get the correct thermal relic density. 

Collider Signatures. The S'-mesons would be pro- 
duced at LHC mainly through intermediate s-channel 
photons and Z bosons in the case where their color per- 
tains to an exotic SU(N) gauge symmetry, as shown in 
fig. 5. Because of confinement the initially produced S-S* 
pair must hadronizc to form the S'-meson bound states 
4> s (SS*) or n s (ST*). If m T < m s (but not <C m s ), 
then 4>s and tts pairs will form with roughly equal prob- 
ability. In the case of tts production, if the S lifetime is 
short enough so that it decays within the detector, there 
will be a distinctive signal of pairs of charged leptons 
and antilcptons, where each like-sign pair has an invari- 
ant mass equal to that of the tts- 

If on the other hand a <j>s pair is produced, the typical 
decay products will be dark matter pairs, XX, or glue- 
balls, ClCl, as shown in fig. 5. Each glucball decays into 
two photons because of the one-loop Eulcr-Hcisenbcrg in- 
teraction. Thus pairs of photons with invariant mass itiq 
will be produced, and pairs of pairs will reconstruct to 
invariant mass m l p s = 2ms- In the example shown, this 
will be accompanied by missing energy m r j )s A detailed 
study should be done to see if existing searches of this 
nature [28] already exclude our model. Another likely 
signature would be four pairs of photons due to decay of 
both mesons into glueballs. 

Conclusions. We have shown that scalar dark mat- 
ter X with mass 130 GeV could produce a gamma ray 
spectral feature tentatively identified in Fermi-LAT data, 
with the addition of just one scalar multiplet S trans- 
forming as (Y s , 1,3) under U(l)y xSU(2) L xSU(3), where 
the SU(3) might be the gauge group of QCD, or else some 
new hidden sector interaction. The coupling XxX \S\ 
between S and X need not be very large unless qs ~ 1 
and ms/rnx is greater than a few. The strong interac- 
tions of S serve two purposes: they confine the stable 
charged relic component of S, and the number of col- 
ors helps to increase the XX — > 77 cross section while 
keeping Xx reasonably small. 

Because S has similar quantum numbers to right- 
handed squarks, it is tempting to make this identifica- 
tion, but such light squarks are ruled out by LHC with 
the possible exception of the third generation. Squarks 
have a smaller electric charge than S in our preferred 
examples, and all three generations would need to con- 
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FIG. 5: Production and decay of S'-mesons at hadron collider. 

tribute to compensate for the resulting decrease in the 
XX —> 77 cross section. 

It is interesting that we rely upon the h 2 X 2 coupling 
between X and the Higgs boson to get the thermal relic 
density of dark matter, and that the same coupling leads 
to a cross section for X scattering on nucleons that is just 
a factor of 3 below the 2012 XenonlOO direct detection 
limit. The model has additional links to Higgs physics: 
the possibility of increasing the h — > 77 branching ratio 
by a factor of a few, and the existence of bound states 
of S and S* (<fis), which could have a large branching 
ratio to decay into Higgs bosons, though more gener- 
ically their decay products are dominated by glueballs 

of the exotic SU(N) or X bosons. We suggest that 
LHC might discover the S'-mesons, whose mass should 
be > 260 GeV, by observation of two photon pairs each 
with invariant mass of mn, accompanied by the same 
amount of missing energy, or four photon pairs each of 



mass mr>. In addition, lighter charged ST* mesons (tts) 
should be produced, decaying into like-sign lepton pairs, 
which might or might not occur within the detector. A 
more detailed study of LHC signals is contemplated [30] . 

After releasing version 1 of this work, we were informed 
of a similar model in [29], resembling ours in the case 
where mx = ni^, s , using the process XX — > 4>s4>S fol- 
lowed by 4>s ~ * 77: which becomes the dominant decay 
channel for <j)g if <C 1 and no hadronic channels are 
available (as in the case of a hidden SU(N) with glueball 
mass greater than m$ s /2). This seems to be another 
viable region of parameter space for our model, relying 
only upon tree-level amplitudes. In this case the DM 
mass would be 260 rather than 130 GeV. We acknowl- 
edge A. Ibarra for pointing this out. 
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